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Abstract 



We describe the dynamical behavior found in numerical solutions of the Vector 
Complex Ginzburg-Landau equation in parameter values where plane waves are 
stable. Topological defects in the system are responsible for a rich behavior. At low 
coupling between the vector components, a frozen phase is found, whereas a gas- 
like phase appears at higher coupling. The transition is a consequence of a defect 
Q\ \ unbinding phenomena. Entropy functions display a characteristic behavior around 

■ the transition. 

O \ 1 Introduction 
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Spatially extended nonlinear dynamical systems display an amazing variety 
of behavior including pattern formation, self-organization, and spatiotempo- 
ral chaos [1]. Transition phenomena between different kinds of states share 
some characteristics with phase transitions in equilibrium systems. Symmetry 
breaking, topological defects, and Goldstone modes, for instance, are com- 
monly found. Nevertheless, a much larger variety of collective effects are pos- 
sible in these far-from-equilibrium systems. 

In this paper we report some numerical results on the behavior of the Vector 
Complex Ginzburg-Landau (VCGL) equation [2,3], a model originally devel- 
oped in the study of pattern formation in optical systems [4]. It consists of a set 
of two coupled complex Ginzburg-Landau equations which could be thought 
as the two components of a vector equation : 

d t A± = A± + (1 + ia)V 2 A ± - (1 + i(3){\A±\ 2 + ^\A T \ 2 )A ± . (1) 



The VCGL equation appears naturally in situations where a two-component 
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vector field starts to oscillate after undergoing a Hopf bifurcation. This is the 
case of the transverse electric vector field in a resonant optical cavity near the 
onset of laser emission. The two complex fields A± are the complex envelopes of 
the two components (the circularly polarized components in the optical case) 
of the oscillating field. The parameter a measures dispersion or diffraction 
effects whereas (3 is a measure of nonlinear frequency renormalization. 7 is the 
coupling between the components, so that for 7 = one obtains two uncoupled 
scalar Ginzburg-Landau equations. 

The onset of oscillations breaks two continuous symmetries. On the one hand, 
the phase of the oscillations destroys time translation invariance. On the other 
the direction of oscillations breaks isotropy by singling out a vector orientation. 
Typically these symmetries are broken differently in different parts of the 
system, so that regions in different oscillation states, with topological defects 
between them, appear and compete. 

For the case 7 = 0, equivalent to the scalar phase diagram charting the 

different states at different parameter values has been obtained both in one 
and in two dimensions [5]. In the general vectorial or coupled case, however, 
our knowledge is much more partial. Here we will describe states appearing 
in two spatial dimensions for 7 real, < 7 < 1, and a and (3 such that plane 
waves A± = Q±e l ^ k± ' yi ~ ul±t ^ are linearly stable solutions (a necessary condition 
is l + a/3 > 0). The range of parameters that we consider is relevant to describe 
laser emission when atomic properties favor linear polarization in a broad area 
laser with large detuning between atomic and cavity frequencies. 

The following two sections describe our results for the behavior of the system in 
our range of parameter values. We show the existence of a transition between a 
frozen phase and a gas-like phase. After the conclusions section, an Appendix 
gives some details on the numerical algorithm used. 



2 Defect-dominated frozen phase 

Despite the existence and stability of plane-wave solutions, typical evolution 
starting from random initial conditions leads to complex evolving states. For 
7 small the state of each component superficially resembles the one obtained 
for the scalar equation (see Fig. 1): the dominant objects are spiral waves, 
emanating from or sinking into a defect (a zero of the complex field, giving a 
phase singularity) core. Despite the similarities, there are important differences 
between defects in the scalar case and in the present vectorial case. In the scalar 
case there is only one complex field, so that there is a single phase and thus 
a single type of charge associated to its singularities or defects. In our case 
there are two complex fields, A + and A_, which can vanish independently, 
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giving rise to two independent charges. The topological charges of a defect are 
defined by 



where T is a closed path around the defect, and the phases 0± are defined by 
the relations A± = \A±\e l< ^ ± . Numerically we do not find in the system the 
spontaneous emergence of any topological charge with modulus greater than 
1 starting from random initial conditions. 

It is possible to make a classification of defects using standard topological ar- 
guments [6,7]. We call vectorial defect a defect which is a singularity of both 
components of the field (i.e. both components vanish at the same point). A 
vectorial defect is of argument type when the charges of the two field compo- 
nents have the same signs, i.e., when n + — n_ = 1 or n + — n_ = — 1. If the 
charges are of opposite signs, i.e., when n + = — n_ = 1 or n + = —n_ — — 1, 
the vectorial defect is of director type. We call mixed defect a defect that is 
present just in one component of the field. 

For 7 = there is no interaction between the two fields and thus binding 
of mixed defects to form vectorial defects would not occur generically. By in- 
creasing 7 we observe that all kinds of defects appear leading to configurations 
which evolve very slowly in time. Such configurations are representative of a 
frozen or glassy state. For example, for 7 = 0.1, a = 0.2, (3 — 2 and large times 
(Fig. 1), the system evolves into a state in which the fields are organized in 
domains of nearly constant modulus separated by shocks. There is a vectorial 
defect at the center of each domain. This defect core emits or receives phase 
waves which entrain the whole domain. Perturbations and mixed defects are 
ejected away from the defect core with a group velocity. The mixed defects 
accumulate at the domain borders. In Fig. 1 we also show the global and 
relative phases, 4> g = <p + + 4>- and <f> r = <p + — </>_. An argument defect has 
a global phase (j) g that rotates An around the defect core, while the relative 
phase <p r rotates 0. For a director defect, 4> g rotates and (f> r rotates 4-k. In 
consequence argument and director defects are easily distinguished in the plot 
of the global phase: A two- armed spiral is formed around an argument defect, 
while a target pattern is seen in the domain of a director defect. Mixed defects 
appear as points around which the global or relative phase rotates by 2n. The 
modulus of this kind of configuration evolves very slowly in time, so that we 
could call it a frozen or glassy state. 

As 7 increases the structure of the mixed defects becomes such that a maxi- 
mum in the modulus of one of the components appears where the other com- 
ponent presents a singularity (see Fig. 2). Such anticorrelation, which also 
occurs for the shocks separating the regions dominated by a vectorial defect, 
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Fig. 1. Frozen configurations for 7 = 0.1, a = 0.2 and (3 = 2, displaying the different 
kinds of defects, (a) |^4+| 2 , (b) |^4-| 2 , (c) global phase <p g , and (d) relative phase 4> r . 
In the modulus plots, black points are zeros of the corresponding field and white 
points locate the maximum values 



becomes more evident by further increasing 7. This feature is also present in 
the one-dimensional case [8]: no topological defects exist in d — 1, but a spa- 
tially localized minimum of one field, which moves in time, goes together with 
a maximum of the other field. 
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Fig. 2. Field configurations for a = 0.2, and (3 = 2. Gray code as in Fig. 1. Top: 
7 = 0.2; wave domains dominated by vectorial defect cores are still present. Maxima 
in the modulus of one component are clearly associated with mixed defects in the 
other. Bottom: 7 = 0.8; defect unbinding has occurred and only mixed defects, with 
associated maxima in the non-singular component, are present. First column: |^4+| 2 , 
second column: |^4-| 2 . 

3 Unbinding transition to a gas phase 

There is a critical value (7 ps 0.35 for a = 0.2, (3 = 2, as in Figs. 1 and 2) 
above which vectorial defects disappear. We observe two different annihila- 
tion processes: a) One of the two singularities that form the vectorial defect 
is annihilated in the collision with a mixed defect, in the same component 
but of opposite charge, which migrates from the boundaries of a domain. A 
mixed defect, with charge associated to the other component, is thus left in 
the system, b) The vectorial defect splits into two spatially separated mixed 
defects, one in each component. 

When the vectorial defects disappear, spiral-wave domains dissolve and the 
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frozen structure transforms into a mobile configuration with fast active dy- 
namics. Fig. 2 (bottom) shows a typical snapshot: mixed defects travel freely 
around the system as in a kind of "gas phase". The anticorrelation between 
the two components is quite evident at this large value of 7. The transition 
between the frozen and the gas behavior is rather sharp, and can be thought 
as a kind of vortex unbinding. 

One way of characterizing the different kinds of behavior and transitions be- 
tween them is by means of an entropy measure H(X) = —J2 x p( x )^ n P( x )y 
where p(x) is the probability that X takes the value x. H(X) measures the 
randomness of a discrete variable X. We can compute the single-point en- 
tropies of the modulus of the field components by considering the discretized 
values of \A + \ and as random variables (X = \A + \ or \A-\; we discretize 
the range of these variables into 200 values). The associated probability dis- 
tributions are defined from the ensemble of values collected from different 
space-time points. In Fig. 3 we plot the entropy of \A + \ and \A_\ as functions 
of 7. For low values of 7 the system is in the frozen state consisting in large 
domains of uniform modulus surrounding vectorial defects. These domains im- 
pose some degree of order which gives low values to the entropies. For 7 = 0.25 
the size of the domains diminishes, and the system becomes more disordered 
as indicated by the increase of the entropies. There is a maximum of the en- 
tropies at 7 ~ 0.3, which is the value at which the argument defects are seen 
to annihilate. Thus the maximum in the entropies is signaling the transition 
from the frozen structure to the gas-like phase. For 7 ~ 0.35 the director vec- 
torial defects disappear also, so that for higher values of 7 there are only mixed 
defects. When 7 leaves the transition region, the entropies initially decrease, 
but they increase later with growing 7 in correspondence with the increasing 
dynamic disorder in the fields. This behavior of the entropies is in contrast 
with the one-dimensional case [8], where topological defects are absent. There, 
entropies maintain an essentially constant value when 7 varies. The presence 
of defects in the two-dimensional case is responsible for the distinct behavior 
of the entropies. 



4 Conclusions 



In this Paper we have described qualitatively some aspects of the dynamics 
of the VCGL equation, focusing in a particular parameter regime of relevance 
in optics. The presence of different kinds of defects is the characteristic phe- 
nomenon organizing other features of the dynamics. Two main "phases", a 
frozen or glassy state and a more dynamic gas-like phase, have been identi- 
fied. The transition between these two phases originates in the unbinding of 
vectorial defects. 
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Fig. 3. Entropy of \A + \ (squares) and \A-\ (triangles) as functions of 7 (a 
(3 = 2). 
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A Numerical Integration Scheme 



The time evolution of the complex fields A±(x, t) subjected to periodic bound- 
ary conditions is obtained numerically from the integration of the VCGL in 
Fourier space. The method is pseudospectral and second-order accurate in 
time. It is the straightforward generalization to two dimensions and two com- 
ponents of the algorithm described in [9] for the scalar Ginzburg-Landau equa- 
tion. Each Fourier mode A q ± evolves according to: 



d t Al{t) = -a q Al(t) + , 



(A.l) 



7 



where a q is (1 + ic\)q 2 — 1, and $± are the g-modes of the non-linear terms in 
the VCGL equation. 

When a large number of modes q is used, the linear time scales a q can take a 
wide range of values. A way of circumventing this stiffness problem is to treat 
exactly the linear terms by using the formal solution: 

A q ± (t) = e- a « l ^(t )e a ^° + J $^(s)e a ^sj . (A.2) 
From here the following relationship can be obtained: 

i _ -2a q 8t 

A q ± {n + 1) = e- 2a " 5t A q ± {n - 1) + - — $Un) + 0{5t 3 ) . (A.3) 

a q 

Expressions of the type f{n) are shortcuts for f(t = n5t). Scheme (A.3) alone 
is unstable for the VCGL equation. To fix this one can derive the auxiliary 
expression 

A q ± {n) = e- a « 5t A q ± {n - 1) + $£(n - 1) + 0(5t 2 ) , (A.4) 

and the algorithm proceeds as follows: 

(1) Starting from A±(n — 1) and Fourier inverting to get A±(x, n — 1) one can 
calculate the nonlinear terms in direct space and then obtain $j_(ri — 1). 

(2) Eq. (A.4) is used to obtain an approximation to A±(n). 

(3) The non-linear terms $±(n) are now calculated from these A±(n) by 
going to real space as before. 

(4) The fields at step n + 1 are calculated from (A.3) by using A±(n — 1) and 
& ± (n). 

At each iteration, we get A q (n + 1) from A q (n — 1), and the time advances by 
2St. 

The number of Fourier modes depends on the space discretization. We have 
used dx — 1 in lattices of size 128 x 128 or 256 x 256. The time step was 
usually dt = 25t = 0.05. 
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